Submitted to Physics of Plasmas 

Analytical description of stochastic field-line wandering in 

magnetic turbulence 

A. Shalchi & I. Kourakis 
Institut fur Theoretische Physik, Lehrstuhl IV: Weltraum- und Astrophysik, 
Ruhr-Universitdt Bochum, Bochum, Germany 

(Dated: February 5, 2008) 

Abstract 

A non-perturbative nonlinear statistical approach is presented to describe turbulent magnetic 
systems embedded in a uniform mean magnetic field. A general formula in the form of an ordinary 
differential equation for magnetic field-line wandering (random walk) is derived. By considering the 
solution of this equation for different limits several new results are obtained. As an example, it is 
demonstrated that the stochastic wandering of magnetic field-lines in a two-component turbulence 
model leads to super diffusive transport, contrary to an existing diffusive picture. The validity 
of quasilinear theory for field-line wandering is discussed, with respect to different turbulence 
geometry models, and previous diffusive results are shown to be deduced in appropriate limits. 

PACS numbers: 47.27.tb, 96.50.Ci, 96.50.Bh 
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I. INTRODUCTION 



Understanding turbulence occupies a central part of current research efforts in space 
physics and astrophysics; see, e.g., in Refs. . Charged particles in turbulent magnetized 
environments perform a complex motion, which may be seen as the sum of a deterministic 
trajectory (helical motion around the magnetic field lines) and a random component, due 
to turbulence. In an effort to elucidate transport mechanisms in turbulent (collisionless) 
magnetized plasmas, some progress has been marked by employing a statistical description 
of turbulence (e.g. [3|, U, [5j ) . The erratic component of charged particle motion can physically 
)e associated with the stochastic wandering (random walk) of magnetic field-lines; see Refs. 



aflflfly,n r 



112, 



13l | . It is nevertheless understood that a definite, widely applicable 



analytical formulation of FLRW (field-line random walk) is not yet available. 

A quasi-linear approach for field-line random-walk formed the model basis in early works; 
see ,e.g., Ref. {14J]. In that description, the unperturbed field-lines are used to describe field- 
line wandering by using a perturbation method. This approach is believed to be correct 
in the limit of weak turbulence, where turbulent fields are assumed to be much weaker 
than the (uniform) mean field (SBi <C B ). A non-perturbative statistical description of 
field-line wandering was later suggested in Ref. 15), relying on certain assumptions about 
the properties of the field-lines (e.g. Gaussian statistics) in combination with an explicit 
diffusive hypothesis for the field-line topology. 

The main scope of this article is to address the problem of field-line wandering analytically, 
in relation with general turbulence models. The limits of the validity of quasilinear theory 
(QLT) are also to be discussed. The standard statistical description of field-line wandering 
is rigorously shown to lead to an ordinary differential equation (ODE) for the mean square 
deviation (MSD) of the field-lines. The ODE is solved in certain cases and the results are 
compared with QLT and associated methods. An interesting example is field-line random 
walk in two-component turbulence, where a super diffusive behavior of field-line wandering 
is clearly found. 
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II. AN ODE FOR FIELD-LINE WANDERING 



We shall consider a collisionless magnetized plasma system which is embedded in a uni- 
form mean field (B Q = B e z ) in addition to a turbulent magnetic field component in the 
transverse direction (SB Z = 0). The field-line equation in this system reads dx/dz = 5B X /B Q . 
Following the established Kubo statistical formalism for random processes, the field-line (FL) 
mean square displacement (MSD) can be written as 

{(Ax(z)) 2 ) = -L f dz' f dz" R xx (z',z") (1) 
-°o Jo Jo 

where we have employed Ax = x(z) — x(0) and the xx— component of the magnetic correla- 
tion tensor R xx (z',z") = (SB x (x(z'))SB*(x(z")))\ the real part of the right-hand side (rhs) 
will be understood. Assuming homogeneous turbulence R xx {z',z") = R xx (\z' — z"\), one 
obtains 

<(Ax(^)) 2 ) = 4 f dz' (z - z')R xx {z') (2) 
n o Jo 

with Rx X (z') = (5B x (x(z'))6B*(x(0))). Since the correlation tensor is itself dependent of 
x(z) = (x(z),y(z), z) this is an implicitly nonlinear integral equation for the magnetic field- 
line space topology. In order to evaluate Eq. ([2]), one has to determine the correlation tensor 
R xx . A Fourier transformation leads to 

R^(z) = Jd 3 k J d 3 k' (^5B x (k)5B:(ky^- ig < ^ , (3) 

or, adopting Corrsin's independence hypothesis [16] 



Rxx(z) = j d 3 k j d 3 k (8B x {k)5B* x [k 

x / e ik-x(z)-ik -x(0)\ _ 



Assuming that the magnetic fields for different wave vectors are uncorrelated 
5B x (k)5B*(k')) = P xx (k)5(k - k') one is led to 



Rxx(z) = j d'kP xx (k){e lk ^) (5) 

with P xx (k) =< 5B x (k)5B*(k) >. For the sake of analytical tractability, in order to evaluate 
the characteristic function < e lk ' Ax ^ >, we shall assume Gaussian statistics for the field- 
lines, thus 

^ e ik-Ax(z)^ = e -i((Ax( 2 )) 2 )fc2-i((Ay( Z )) 2 )fc2+ifcn^ ^ 



For axisymmetric turbulence ((Ax) 2 ) = ((Ay) 2 ), so Eq. (JSJ) takes the form 

((Ax(z)) 2 ) = / d 3 kP xx (k) 



x / dz' (z-z')cos(k\\z')e-^ Ax{z '^) k ^. (7) 
Jo 

Upon differentiation with respect to z, we find for the field-line MSD 

^<(A^)) 2 > = Jj jTW ^00 
<i 3 fc P xx (k) 

x / ^ cos(A;||Z / )e"^ (A ^ ,))2 ) /ci . (8) 





A second differentiation leads to 

<Z 2 



2 f 3 
= ~~j^2 J d k P xx (k) 



i((Ax( 2 )) 2 )fei_ ^ 



x cos{k\\z)e 2 

This ODE was obtained, relying on no other assumptions than Corrsin's hypothesis and 
Gaussian FL statistics. In combination with different turbulence models, it provides a gen- 
eral basis for the determination of the FL-MSD, thus allowing for a quantitative description 
of field-line wandering. In the following, we shall consider a specific example, as well as 
various limits of this description. 

III. ANALYTICAL AND NUMERICAL RESULTS FOR SLAB/2D COMPOSITE 
TURBULENCE 

A two-component turbulence model has been proposed as a realistic model for solar wind 



turbulence [l7j. Within this model, the turbulent fields are described as a superposition of 
a slab model (k \\ Bq) and a two-dimensional (2D) model (k _L -Bo)- In the following, we 
shall evaluate the field-line MSD, in view of identifying an comparing among pure slab and 
composite slab/2D geometry. 
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A. Analytical results for pure slab geometry 



A first approach consists in assuming slab turbulence statistics. The 
xx— component of the correlation tensor then reads 



P^(k)=g^\k ll ) 6 -^, (10) 



where we assume 



9 sla \k\\) = ^r-lsi ab 6B 2 slab (l + k\ll ab r (11) 
for the slab wave spectrum. Here C{y) is a normalization constant given by 

^sW^b) (12) 

(see e.g. [l8l|). l s i a b is the slab-bendover-scale, 2v is the inertial-range spectral 
index, and 5B 2 lab / B^ determines the relative strength of slab turbulence. For the 
spectrum of Eq. ( 1111) . it is straightforward to show that the slab result behaves 
diffusively, as 

— ^^slab | % | (13) 

with the slab field-line diffusion coefficient 



8B 2 

Kslab = 7cC(v)l slab —g±. (14) 

This is an exact result, readily obtained upon analytical evaluation of Eq. ((§]) 
in the limit I z I— > oo. 



B. Analytical results for slab/2D composite geometry 

In a hybrid (composite) slab/ '2D model, the correlation tensor is assumed to 
have the form P xx (k) = P^ b {k) + P 2 ®{k) (following the definitions above). Here 
we used the slab correlation tensor of Eq. ( 1101) and the 2D correlation tensor 
defined by 

^) = ^W^(l-§). (15) 

Eq. ((9j) thus becomes 

+ % rdk ± g™(k ± )e-H^) 2 K. (16) 
-°o Jo 
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Focusing on high values of the position variable | z \— > oo, hence ((Ax) 2 ) — > oo, the 
main contribution to the integral comes from very low values of the integration 
variable. Assuming a quasi-constant behavior of the 2D spectrum in the energy- 
range we may approximate as 

+ %g 2D (V) rdk ± e-^m. (17) 
-°o Jo 

It can easily be demonstrated, by evaluating Eq. ([3]) for pure slab geometry, 
that the slab correlation function R s ^ b for the spectrum of Eq. ( Illj) has the 
asymptotic behaviour 

ie 6 (M-oo)~ feV Vw/'-». (18) 



, z 

Obviously the first contribution in Eq. ( 1171) is much smaller than the second, 
thus one obtains 

dz B V 2 ^(Ax) 2 ) 

It is straightforward to show that Eq. ( 1191) is solved by 

-i 2/3 



;(Ax) j > 



(20) 



in the limit | z \— > oo. To proceed we use a form for g 2D (k±) similar to the slab 
spectrum (see Eq. (1111 )) 

g2D = 2C(v)_ 12d5B 2 (1 + tfp y r (21) 

7T 

Here 1 2 d is the 2D-bendover-scale, 2u is the inertial-range spectral index, and 
bB\ D jBl determines the relative strength of 2D turbulence. Combining this 
model spectrum with the above, we find 

2/3 



((Ax) 2 ) 



9C(v)J-l 2D 



7T. 6Bi D 



2- Bl j 



z | 4 / 3 (22) 



which is clearly a non-diffusive result. As demonstrated, field-line wandering 
behaves super diffusively if the slab/2D composite model is employed for the 
turbulence geometry. 



C. Numerical evaluation 



Eq. (J7|) is an integral-equation for the MSD of the field-lines. By using Eqs. (|T5|) and 
(|2"T|) for P 2 ®(k), in addition to Eqs. (fTOj) and ffTTl) for P* l x ab (k), the total correlation tensor 
Pxxik) for the composite model is specified. 

We have evaluated Eq. ([7]) numerically, adopting a 20% slab-/80% 2D- composite turbu- 
lence model. For the sake of comparison, we have also computed the result in the pure-slab 
turbulence case for 5B 2 lab /B 2 = 1. For the same cases of study, we have seen above that the 
analytical result for two-component turbulence is given by Eq. (1221) . whereas the pure slab 
result is given by (TT^|) . 

The numerical results are compared with the analytical results (presented above) in Fig. 
[U by depicting the running diffusion coefficients < (Ax(z)) 2 > /(2zl s i a b) as a function of 
z/hiab- An excellent agreement is witnessed among the theoretical and the numerical results, 
in both cases. The analytical expression (122]) derived above for the asymptotic behavior of 
the field-line MSD thus appears to be valid in the composite turbulence model also and, 
in fact, bears a contributions which soon exceeds the diffusive result of the slab-model 
significantly: cf. the upper curve(s) in Fig. 1 to the lower, constant curve(s). 

IV. FURTHER RESULTS AND LIMITS 

We shall complete our study by considering various limiting cases, in which previous 
theoretical results are recovered. Although these results are not new, it is interesting per se 
to demonstrate that they can be obtained as appropriate limits from Eqs. ([319]). 

A. The initial free streaming regime 



For small values of the position variable (z — > 0) hence expecting ((Ax) 2 ) — > 0, we obtain 
from Eq. © 



Assuming vanishing conditions for both the MSD and its derivative at zero, one obtains 




(23) 




(24) 
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so a strong superdiffusive FL wandering regime is clearly found for small z, independent of 
the turbulence model adopted. For pure slab geometry it can easily be demonstrated 
that the initial free streaming solution is valid for | z |<C I slab- This means that 
Eq. ( 1241) is valid for length scales which correspond to the inertial range of the 
turbulence wave spectrum. 



B. Quasilinear theory for field-line random walk 

Within QLT, one may replace the MSD on the right hand side of Eq. by the unper- 
turbed field-lines (thus taking ((Ax) 2 ) = 0), which yields 

^ ((Ax) 2 ) = ^iJd 3 k P xx {k) coB(k\\z). (25) 

QLT is thus apparently only exact for pure-slab turbulence, where P\ m ~ S(k±); cf. Eq. 
( TTUj) . On the other hand, within QLT one finds for pure 2D and slab/2D composite tur- 
bulence [adopting Eqs. ( |T5l) and ( 125]) ] that ((Ax) 2 ) = (SB X / B ) 2 z 2 , in disagreement with 
the nonlinear result obtained above. Thus, QLT fails to describe field-line wandering in the 
two-component model. Presumably, quasilinear theory might thus also not apply in other 
non-slab models. 



C. The diffusion limit 

Combining Eq. (jSj) with the assumption of diffusion for the magnetic field-lines (i.e. 
explicitly setting ((Ax) 2 ) = 2k | z |) one obtains 

K = \ [ d 3 k P xx (k) [ dz cos(h { z)e- KZk ±. (26) 
5 o J Jo 

The z— integral can easily be solved to give 

K ^/ ft «>i[W' (27) 

This formula is correct if field-line random walk behaves diffusively and if the small length 
scales of the initial free streaming regime are unimportant. 

For two-component turbulence Eq. (|2"7|) can easily be evaluated; we find 



k 2 

K = K slab + , (28) 

K 



using the slab diffusion coefficient n s iab 

of Eq. (P3J and 

k1d = ^iI d " k kl2p ™ {%) • (29) 

Eq. fl28l) is a quadratic equation in k, which may easily be solved to get 

K _ K slab + \Z K Lfc + jjjjfp ^gQ^ 



We note that this coincides with the result derived earlier by Matthaeus et al. 15]. For 
pure slab turbulence, we have k 2 d = 0, so the expected limit n = K s i a b is recovered. On the 
other hand, for pure 2D geometry, i.e. for n s i a b = 0, one finds k = k,2D- Thus the parameter 
k,2d can be identified with the "diffusion coefficient" for pure 2D turbulence. However, 
adopting the standard spectrum of Eq. ( 1211) . one clearly obtains a diverging 
result: 

r f)2 poo 

k\ d = 2C(v)l 2D d -^f / dk ± k-* {I + lj D kiy u 
£>o Jo 

-> oo. (31) 

The only possibility to prevent this singularity is the introduction of a finite 
box-size L 2 d 3> Iid of the 2D fluctuations. Then we find 



"2D — 2C(z/)/ 2 D 



dk ±k f(i+n D kir 

^0 JL-Zh 



J 2D 

r, ( 1 3 L2 \ 

x 2 F 1 f + v + -;-_) ' ( } 

where the integral was solved by applying exact relations which can be found 



in Ref. |191 ]. Because of L 2 d 3> ho we consider the asymptotic limit of the 



hypergeometric function 2 Fx(a,b; x) (see [20]) to find for k 2 d the expression 

K2D = 



^ VI2DL2D p 2D - (33) 



2v + 1 v B 

Within the diffusion theory for FLRW, the diffusion coefficient is strongly con- 
trolled by the box-size L 2 d- In the light of the results presented in section III., 
the true reason for the divergent behavior (see Eq. ( 131ft ) may be the superdif- 
fusive nature of field-line random walk. It appears that Eqs. ( 1301) and ( 133ft fail 
to provide an appropriate description of the composite turbulence, in contrast 
with the generalized nonlinear theory presented in this paper. 



V. IS A SUPERDIFFUSIVE BEHAVIOUR OF FLRW REASONABLE? 



A key result of the current article is the super diffusive behaviour of FLRW for 
the slab/2D composite turbulence model (see Eq. ( 1221) and Fig. 1). However, 
this nondiffusive result is clearly in disagreement with the assumption of diffusive 
field line wandering employed in several previous articles (see e.g. {15I]). 

The super diffusive result derived here is based on two ad hoc assumptions 
which cannot be deduced from first principles: 

1. A random phase (or Corrsin type) approximation was applied (see Eq. ([4])). 

2. A Gaussian distribution of field lines was assumed (see Eq. (IHD). 

Eq. ( 1221) therefore provides an accurate result, suggesting a super diffusive be- 
havior of FLRW in real systems, provided that the latter two assumptions hold. 
However, one might draw the conclusion that the non-diffusivity deduced in this 
article is a consequence of (two) inaccurate approximations. We shall attempt 
to remove this ambiguity in the discussion that follows. 



In another article ([2l|]) we show that already for pure slab fluctuations where 
FLRW can be described without any assumptions or approximations, a superdif- 
fusive result can be obtained be replacing the standard spectrum of Eq. ( 1111) 
by a decreasing spectrum in the energy range. According to these results, su- 
perdiffusion indeed arises as the normal, by default behavior of FLRW. Classical 
(Markovian) diffusion can only be obtained in very special limits (wave spectrum 
exactly constant in the energy range and pure slab geometry). 

The slab correlation function decays rapidly with increasing distance (see Eq. 
( 1181) ). For decreasing wave spectra in the energy range and for non-slab models, 
however, we find a much weaker (non-exponential) decrease of the magnetic 
correlation functions. This weaker descrease of magnetic correlation functions 
provides a physical explanation of the nondiffusivity discovered in the current 
article. 

A useful application but also a test of the super diffusive result is cosmic ray 



scattering perpendicular to the mean magnetic field. In jl8l ] the super diffusive 



result obtained here is combined with a compound transport model based on 
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the assumption that the guiding centers of the charged cosmic rays follow the 
magnetic field lines. In this article a quite good agreement between test particle 
simulations and the combination of compound transport and superdiffusion of 
FLRW was found. This agreement can also be seen as a proof of the validity 
of the two assumption (Corrsin's hypothesis, Gaussian statistics) used in this 
article to deduce the superdiffusive result. 

Furthermore, it is argued in [l8| that for a diffusive behaviour of field line 
wandering one always obtains a strong subdiffusive behaviour of cosmic ray 
perpendicular scattering which cannot be true. Thus, the superdiffusive result 
deduced in the current article is not only a result obtained under certain assump- 
tions. It seems that the real field lines indeed behave super diffusively and that 
this superdiffusivity is essential for understanding charged particle propation. 

From a purely fundamental point of view, the statistical-mechanical descrip- 
tion of FLRW in magnetized plasmas bears certain generic physical characteris- 
tics which are encountered in various physical contexts. Our analytical findings 
are qualitatively reminiscent of earlier results by Isichenko on advection-diffusion 



problems 22]. Physically speaking, diffusive behavior of turbulence implies a 
fast decay of Lagrangian correlations, sufficiently fast for the limit < (Ax) 2 > /2t 
to converge. In fluid (non-plasma) environments, long scale random velocity 
fluctuations (depending on boundary conditions) may dominate the motion of 
a fluid element and may thus cause the limit to diverge. On the other hand, 
turbulent magnetized plasma systems are characterized by a plethora of non- 
linear mechanisms, intertwined with the complex FL topology, which may add 
up in such a way that anomalous (non-classical) diffusion become the rule in 
plasmas. Persisting correlations among successive FL displacements are thus 
built-up, resulting in the superdiffusive behavior exposed in our work. This 
seems to suggest a general tendency of turbulent magnetized systems, so that 
previous diffusive FLRW considerations (obtained under various assumptions) 
should rather be treated with some caution. 

It may be added that a superdiffusive result for FLRW was also obtained in 
various previous works, where distinct methodologies were employed. In Ref. 

, the FL separation across the field was studied by adopting an ad hoc form for 
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the magnetic correlation tensor; its behavior was found to be very sensitive to the 
exact form of the power spectrum, while a parabolic (ballistic) regime was shown 
to rule for steep spectra. In Ref. |7||, geometric (non-statistical) arguments were 
employed to study FLRW in galactic cosmic rays. A super diffusive behavior was 
also found, independently from the turbulence spectrum (though mainly shallow 
spectra were considered therein). Similar findings were reported in Ref. [8j, in 
the context of thermal conduction in galactic clusters. A chaotic magnetic FL 
behavior (explicitly assumed to obey a Lyapunov behavior over a wide range 
of length scales), i.e. assuming an exponential FL divergence, was there shown 
to enhance perpendicular thermal conduction significantly, essentially increasing 
thermal conductivity in the transverse direction (previously thought to be much 
weaker than the parallel one) by a factor (up to) 5, as compared to the Spitzer 
23]. Still in the galactic thermal conduction context, Refs. |1(3, 11] have 



value 



adopted a double diffusion concept, where particle diffusion along the field was 
related to FLRW in space, in order to study thermal conduction and electron 
diffusion along the magnetic field; a combined phenomenological, Fokker-Planck 
and numerical approach, FL displacement was shown to increase fast over length 
scales which are a multiple of the characteristic turbulence scale. Also note Ref. 
Q] (and the exhaustive discussion therein) for a nonlinear treatment of the 
perpendicular diffusion of charged particles (not FLRW). 

VI. CONCLUSION 

A generalized nonlinear formulation has been presented, to describe field-line wandering 
(random walk) in general magnetostatic systems consisting of a statistical (or turbulent) 
component 5B and a nonstatistical uniform component Bq. Assuming vanishing parallel 
component of the turbulent field (5B Z = 0), applying the Corrsin approximation, and as- 
suming Gaussian statistics, a general ODE was deduced for the field-line MSD (Eq. Q)) in 
axisymmetric and homogeneous turbulent plasmas. 

Adopting the two-component turbulence model which was suggested by Bieber et al. 



171 ]) as a realistic model for solar wind turbulence, it was demonstrated systematically 



that FLRW behaves superdiffusively. Specifically, a weakly superdiffusive behavior of the 
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mean square deviation was found in the form < (Ax) 2 >~| z | 4 / 3 . This result was compared 
to previous results obtained via different assumptions, namely the quasilinear result (< 



(Ax) 2 >~ z 2 ) and the diffusive result < (Ax) 2 >~| z | of Ref. [15]. 

Extending the considerations in this article, the generalized nonlinear formulation pre- 
sented here may be applied in a description of non-axisymmetric turbulence and/or other 
(e.g., anisotropic) turbulence geometries. 

Since field-line random walk is of major importance in the description transport of charged 
particles in astrophysical plasmas, these results are significant in cosmic ray space physics 
and astrophysics. Concluding, the analytical expression(s) derived in this article provide a 
useful toolbox, which extends existing transport theoretical models and might be essential 
in inspiring forthcoming ones. 
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FIG. 1: (Color online) The running diffusion coefficients (^(Ax(z)) 2 J /(2zl s i a i,) as a function of 
z /hiab are depicted. The analytical result (solid line) is compared to the numerical result (dashed 
line) for 20 % slab / 80 % 2D composite geometry. The analytical slab result (dash-dotted line) 
and the numerical slab result (dotted line) are also provided, for reference. The analytical and 
numerical results are are obviously in good agreement. 
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